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We study the dynamics of vortex lattice formation of a rotating trapped 
Bose-Einstein condensate by numerically solving the two-dimensional Gross- 
Pitaevskii equation, and find that the condensate undergoes elliptic deforma- 
tion, followed by unstable surface-mode excitations before forming a quan- 
tized vortex lattice. The dependence of the number of vortices on the rotation 
frequency is obtained. 

PACS numbers: 03.75.Fi, 67.40.Db. 
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Quantized vortex lattices in the atomic Bose-Einstein Condensates 
KA , (BECs) have recently been created experimentally by the groups at ENSlil 

and MITB While several paperaJ have discussed equiliblium configurations 
of vortex lattice, this paper describes the nonequiliblium dynamics towards 
them by numerically solving the Gross-Pitaevskii equation (GPE) that gov- 
erns the time evolution of the condensate order parameter ^(r, t): 

(i_^)n^^:M=[_|!.V2 + 14rap(r)+5lV'(r,t)p-/i-l^L,]v^(r,i). (1) 

Here g = Airh a/m is the coupling constant, proportional to the ®^Rb scat- 
tering length a ~5.77nm. The high anisotropy of the cigar-shaped potential 
used in the ENS experiments may permit the two-dimensional analysis. We 
thus focus on the two-dimensional dynamics of Eq.(ffl) by assuming the trap- 
ping potential 

Hrap(r) = i^mujl{{l + e^)x^ + (1 + ey)y^}, (2) 
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where oj±_ = 27r x 219Hz, and the parameters ex = 0.03 and Ey = 0.09 describe 
small deviations of the trap from the axisymmetry, corresponding to the 
ENS experiments.a The centrifugal term —QLz = ihQ,{xdy — ydx) appears 
in a system rotating about the z axis at a frequency 0. An important 
characteristic parameter of the two-dimensional system is C = SirNa/L, 
with the total number N of the condensate atoms and L the size of the system 
along the z axis. The term with 7 in Eq. (^) introduces the dissipation. 

In this paper we investigate the response of the condensate to the sudden 
turn-on of the rotation of the potential. Section 2 describes the scenario of 
the vortex lattice formation, following the elliptic deformation of the conden- 
sate and the surface-mode excitations.a In Sec. 3 we study the dependence 
of the number of vortices on 0, and the critical frequency as a function of 
C. 



2. FORMATION OF VORTEX LATTICE 

We first prepare an equilibrium condensate in a stationary potential. 
Figure 1 shows the typical dynamics of the condensate density |'(/'(r,t)p af- 
ter the potential begins to rotate suddenly with i7 = 0.7u;_l.H The condensate 
is elongated along the x axis because of the small anisotropy of Ftrap, and 
the elliptic cloud oscillates. Then, the boundary surface of the condensate 
becomes unstable, exciting the surface waves which propagate along the sur- 
face. The excitations are likely to occur on the surface whose curvature is 
low, i.e., parallel to the longer axis of the ellipse. The ripples on the surface 
develop into the vortex cores around which superflow circulates. Subject to 
the dissipative vortex dynamics, some vortices enter the condensate, eventu- 
ally forming a vortex lattice. As the vortex lattice is being formed, the axial 
symmetry of the condensate is recovered by transferring angular momentum 
into quantized vortices. 

This peculiar dynamics is understood by investigating the phase of 
V'(r, t) as shown in Fig.^. As soon as the rotation of the trapping poten- 
tial starts, there come in some branch cuts between the phases and 27r, 
whose ends represent phase defects, i.e. vortices. When the defects are on 
the outskirts of the condensate where the amplitude [-(/^(r, i)| is almost neg- 
ligible, they neither contribute to the energy nor the angular momentum of 
the system. These defects come into the boundary surface of the conden- 
sate within which the density grows up and the phase defects compete with 
each other and induce the above surface waves due to interference. Then 
the selection of the defects starts, because their further invasion costs the 
energy and the angular momentum. As is well known in the study of ro- 
tating superfluid helium,l3 the rotating drive pulls vortices into the rotation 
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Fig. 1. Development of the condensate density |^p after the trapping po- 
tential begins to rotate suddenly with fi = 0.7uj± for C = 1500. The time is 
t =0 msec(a), 21 msec(b), 107 msec(c), 114 msec(d), 123 msec(e), and 262 
msec(f). The unit for length is ano = \/h/2mujj_ = 0.512/im. 
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Fig. 2. Phase profile of ip. The value of the phase varies from (black) to 
27r (white). The time of each figure is the same as that in Fig|l| The unit for 
length is the same as that of Figjl|. 
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Fig. 3. Time evolution of the distortion parameter a (solid line) and the 
angular momentum izl'h per atom (dashed line) in the dynamics of Fig.l. 



axis, while the intervortex repulsive interaction tends to push them apart; 
their competition yields a vortex lattice whose vortex density depends on 
the rotation frequency. In our system, some vortices enter the condensate 
and form a lattice dependent on fi, while excessive vortices are repelled and 
escape to the outside. Remarkably, the phase profile of Fig.2(f) reveals that 
the repelled vortices also form a lattice on the outskirts of the condensate. 
Since they cannot be seen in the corresponding density profile of Fig. 1(f), 
they may be called "ghost" vortices. 

The distortion of the condensate to an elliptic cloud was theoretically 
studied by Recati et al.Q Assuming the quadrupolar velocity field v(r) = 
oSJ{xy\ they obtained the distortion parameter a = ^^(i?^ — i?^)/(i?^ + i?^) 
in the steady states as a function of i7, where Rx and Ry are the sizes of 
the condensate along each direction. Madison et al. observed that, after the 
rotation of 17 = 0.7u;j_ starts suddenly, a/uj±_ oscillates during a few hundred 
milliseconds and then falls abruptly to a value below 0.1 when vortices enter 
the condensate from its boundary surface.Q Figure 3 shows the oscillation of 
a and the increase of the angular momentum ^^/Ti per atom in our dynamics 
of Fig.l. This figure closely resembles Fig. 3 of Ref. 7, and our scenario is 
consistent with the experimental results. 

The dispersion relation of the surface wavesH of a rotating condensate is 
given by uP' = R{uj\ — V!?')k with the condensate size R. The dependence of w 
on R shows that the surface with lower curvature is excited more easily, which 
is clearly seen in Fig.l. The group velocity duj/dk = J R{uj\ — Vt'^)/{2\/k) 
agrees with the propagation velocity of the surface wave in our simulation. 
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3. CRITICAL FREQUENCIES AND MULTIPLE- VORTICES 

LATTICES 

We have examined the critical frequency at which vortices can enter. 
Figure 4(a) shows the dependence of the number of vortices on the fre- 
quency Vt,/uj±, for C =250, 450 and 1500. For C = 450 corresponding to the 
experimental conditionjll we found that only when Q. exceeds Q.ci — 0.56u;j_ 
vortices can enter the condensate and form a lattice. This critical frequency 
is closer to the observed value Onuc = 0.65tj_L, than the values obtained in 
previous literature.Ej The increase in C reduces the critical frequency, and 
stabilizes the lattices of more vortices for the same frequency, The increase 
in n causes lattices of multiple vortices; figure 4(b) shows a lattice of 44 
vortices for C =1500. 

The critical frequencies have been studied from thermodynamic or sta- 
bility arguments ,li30 but they are generally much smaller than the observed 
value. It should be noted that these critical frequencies only give the nec- 
essary condition which enables a vortex to exist stably at the center of the 
trap. Actually, vortices should be nucleated at the boundaries and come into 
the condensate; the condition of $7 which realizes such nonlinear dynamics 
may be generally different from that obtained from the stability arguments. 
Isoshima and Machida examined the local instability of nonvortex states to- 
ward vortex states within the Bogoliubov theory .1111 They note that this local 
instability could give higher critical frequencies than the global instability 
which comes from comparing their total energy. From Fig. 4(a) the critical 
frequencies are 0.58, 0.56 and 0.45, respectively for C =250, 450, and 1500. 
They agree with the critical frequencies 0.57, 0.52 and 0.41 for the same 
values of C obtained from the local stability analysis .1111 



4. CONCLUSIONS 

We have studied the dynamics of vortex lattice formation of a rotating 
trapped BEC by numerically solving the two-dimensional GPE, and obtained 
the following physical picture. When the trapping potential begins to rotate 
at a sufficiently fast frequency, the condensate is distorted to an elliptic 
shape and oscillates. Then its boundary surface becomes unstable, exciting 
surface waves. The origin of these ripples is identified to be the interference 
of violent phase fiuctuations that occur on the outskirts of the condensate; 
and some of the surface ripples develop into the vortex cores which then 
form a vortex lattice. The dynamics is concerned with the "ghost" vortices 
which are the phase defects outside the condensate. 
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Fig. 4. (a) Number of vortices versus n/o;^, for C =250, 450 and 1500. (b) 
Lattice of 44 vortices for C =1500 and Vijuj^ = 0.975. 
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